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qq ■ Abstract 

We apply the new generation of theoretical models of low-mass stars to secondaries in CVs, focussing on 
systems above the period gap. The models confirm that the spectral type should be a good indicator of the donor 
mass. The orbital period-spectral type diagram can potentially constrain the long-term mean mass transfer rate. 
A transfer rate that increases with decreasing period is most easily reconciled with the observational data. 

\Q • 1. Introduction 

o: 

Q\ , The new generation of low-mass star and brown dwarf models by Baraffe et al. (1995, 1997, 1998, 

0^ ' henceforth summarized as BCAH) represent a significant improvement in the quantitative description of 

(— ! \ stars with mass < 1 M . 

The main strengths of the models are in two areas: the microphysics determining the equation of state 
J^ \ (EOS) in the stellar interior, and the realistic non-grey atmosphere models which enter as the outer 

boundary condition. The EOS (Saumon, Chabrier and Van Horn 1995) is specifically calculated for very 
low-mass stars, brown dwarfs and giant planets. It has been successfully tested against recent laser- 
driven shock wave experiments performed at Livermore, which probe the complex regime of pressure 
dissociation and ionization relevant for these objects (cf. Saumon et al. 1998). Recent much improved 
cool atmosphere models (see e.g. the review of Allard et al. 1997) now provide realistic atmosphere 
profiles, which we use as the outer boundary condition, and synthetic spectra. Chabrier & Baraffe (1997) 
have shown that evolutionary models employing a grey atmosphere instead, e.g. the standard Eddington 
approximation, overestimate the effective temperature for a given mass, and yield too large a minimum 
hydrogen burning mass. 

Several observational tests confirm the success of evolutionary models based on these improvements, 
e.g. mass-magnitude relations, colour-magnitude diagrams (Baraffe et al. 1997, 1998), mass-spectral 
type relations (Baraffe & Chabrier 1996), the first cool BD GL 229B (Allard ct al. 1996), properties of 
components in detached eclipsing binaries (Chabrier & Baraffe 1995) and of field M-dwarfs (cf. Beuer- 
mann et al. 1998). 

Here we consider the BCAH models in the context of CV secondaries, and focus on the relation 
between spectral type (SpT), donor mass (M2) and orbital period (P). We obtain the spectral type of a 
stellar model from its calculated colour (/ — K) and the empirical SpT — (I — K) relation established by 
Beuermann et al. (1998). A summary of the input physics used for the calculations presented below is 
given by Kolb & Baraffe (1999). 
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2. The mass spectral type relation 

Mass loss causes stars to deviate from thermal equilibrium. The surface luminosity is no longer in 
balance with the luminosity generated in the core by nuclear reactions, and the difference causes the star 
to expand or contract. Therefore the stellar radius can be either larger or smaller than the corresponding 
equilibrium radius (see e.g. Whyte & Eggleton 1980). Remarkably the effective temperature of low- 
mass main-sequence stars is rather insensitive to the degree of thermal disequilibrium. They behave just 
like giant stars on the Hayashi line, expanding along an evolutionary track with nearly constant effective 
temperature. This can be understood in terms of simple homology relations for predominantly convective 
low-mass stars (e.g Stehle et al. 1996). 

Using BCAH models we quantify the deviation from equilibrium spectral type, as a measure of the 
surface temperature, for CV donors with various evolutionary histories. We consider the following simple 
cases: 

Standard sequence: mass transfer starts from an unevolved (ZAMS) donor with mass 1 M@, proceeds at 
a constant rate 1.5 x 10~ 9 M©yr _1 , stops when the donor becomes fully convective (at mass 0.21 
M Q ), and resumes at the lower rate 5 x 10~ n M Q yr _1 once the donor has settled back into thermal 
equilibrium. This sequence fits the period gap in the framework of disrupted orbital braking (e.g. 
King 1988, Kolb 1996). 

Unevolved sequences with constant M : mass transfer starts from an unevolved (ZAMS) donor with mass 
1 M Q and proceeds at a constant rate M (10 -8 M Q yr _1 or 10 -7 M Q yr _1 ). 

Evolved sequences with constant M : mass transfer initiates from a donor which has already burned a 
significant fraction of its hydrogen supply but is still in the core hydrogen burning phase. We 
show three examples: a moderately evolved low M sequence (initial central hydrogen content 
X c = 0.16, initial donor mass M 2 = 1.0 Mq, M = 1.5 x 10~ 9 M Q yr _1 ), an evolved, low M 
sequence (X c = 4 x 10~ 4 , M 2 = 1.2 M _, M = 1.5 X 10 -9 M^r -1 ), and an evolved high M 
sequence (X c = 4 x 10~ 4 , M 2 = 1.0 M Q , M = 5 x 10~ 8 Moyr" 1 ). 

Note that as the mass loss timescale is short compared with the nuclear time, the nuclear state of the 
donor is essentially frozen once mass transfer has started. 

Fig. [I] shows that the effect of thermal disequilibrium is negligible along the standard sequence; this 
hardly differs from that for CVs with a ZAMS donor. If M is higher than in the standard sequence, the 
spectral type is slightly later for a given mass, while in a sequence with evolved donors it is somewhat 
earlier. As most CVs should form with essentially unevolved donors (Politano 1996, de Kool 1992), and 
as the estimated M exceeds 10 -8 M Q yr _1 only for very few CVs (e.g. Warner 1995, his Fig. 9.8), the 
upper solid curve (X c = 0.16) and the dashed curve [M = 10~ 8 M Q yr~ 1 ) in Fig. [l] should bracket the 
location of the majority of CVs in the mass-spectral type diagram. The figure suggests that the spectral 
type is a good indicator of the donor mass, with a typical uncertainty AM 2 — O.IMq, for any given SpT 
earlier than M0-M2. For later SpT this mass "determination" becomes impractical as the curves in the 
figure flatten and AM 2 is large. Even with extreme evolutionary cases, the unevolved sequence with very 
high transfer rate (M = 10~ 7 M Q yr _1 ) more typical for supersoft X-ray binaries than for CVs, and the 
evolved sequence starting on the terminal main-sequence, the range covered by the secondary stars in the 
SpT-m&ss diagram remains surprisingly narrow (Fig. 0; see Kolb, King & Baraffe 1999 for more details). 

Although the observations compiled by Smith & Dhillon (1998) do not seem to support this prediction, 
the observational errors on M 2 are far too large for any definitive conclusions to be drawn from this 
mismatch. On the contrary with improved and more reliable mass determinations, one may hope to 
constrain the range of M and the relative importance of evolved and unevolved systems in the CV 
population. 
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Fig. 1. Spectral type versus secondary mass. Dots 
denote ZAMS models. Solid curve, close to the dots: 
standard sequence; solid curve, above the dots: donor 
moderately evolved (X c = 0.16). Dashed: donor un- 
evolved, high M (10~ 8 Moyr" 1 ). 



Fig. 2. Same as Fig. |l| but with extreme assumptions 
for the evolved sequence (low M, X c — 4 x 10 -4 ) and 
the unevolved high M sequence (10 -7 Moyr -1 ). 



The figures also apply to secondaries in LMXBs. Evolutionary considerations suggest that, unlike CVs, 
most short-period LMXBs, form with an evolved main-sequence donor with mass > 1 Mq (King & Kolb 
1997, Kalogera et al. 1998). 

3. The orbital period spectral type diagram 

A second invaluable diagnostic diagram for CVs is the orbital period-spectral type diagram (PS di- 
agram) — the "HR diagram" analogue for CV secondary stars. In a semi-detached binary the Roche 
lobe filling star's mean density p determines the orbital period P almost uniquely (e.g. King 1988), 

1/2 

Ph = k/ pg , with k ~ 8.85 being only a weak function of the mass ratio, Ph the period in hr, and p© 
the mean density in solar units. Hence the PS diagram probes the donor structure in terms of two stellar 
parameters that are relatively easy to determine from observation. (The importance of the PS diagram 
has already been pointed out by Ritter 1994). 

The ZAMS track defines an upper bound for the scatter of observed CV secondaries in Fig. Q. Systems 
immediately above and below the period gap have the same spectral type. The results above suggest 
that these CVs have the same donor mass, consistent with the standard period gap model according to 
which CVs evolve through the gap with constant secondary mass. 

The location of various evolutionary tracks above the period gap are shown in Fig. g. Most tracks 
intersect the ZAMS track at a critical period P cr ;t — 5 — 6 hr, which corresponds to a donor with mass 
Mi ~ 0.6 M Q marking the transition from mainly radiative to predominantly convective stars. A simple 
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Fig. 3. Spectral type versus orbital period for CV 
secondaries. Data taken from Beuermann et al. 1998. 
Solid: theoretical location of ZAMS secondaries. 



Fig. 4. Same as Figj3[ but with various evolution- 
ary sequences, cf. also the track morphology in Fig. |s|- 
Thick solid: standard sequence; thin solid: evolved low 
M sequence (X c — 4 x 10~ 4 ); short/long dashes: un- 
evolved sequence (M = lO _8 M0yr _1 ); thick dashes: 
evolved high M sequence (X c — 4 x 10~ ). 



corollary from the results of Sec. 2 is that the donor mass on any evolutionary track is about the same as 
the one on the ZAMS track with the same spectral type. Deviations to shorter or longer periods from the 
ZAMS track mean that the donor is smaller or larger then the corresponding ZAMS star, respectively. 
From the morphology of tracks in Fig. ra it is clear that unevolved high M sequences are excluded at 
long periods (P>P C rit), but would explain the observed scatter for P<P cri t. Likewise, the observations 

exclude evolved low M sequences at short periods, while these are clearly required for P>P cr j t . 

This suggests that, overall, a secular mean mass transfer rate increasing with decreasing orbital period 
would provide a good fit to the observed PS diagram above the period gap. CVs would evolve along 
a low M track at long periods and along a high M track at short periods. To explain the scatter of 
systems at long periods the full spectrum of nuclear evolution prior to mass loss is required. Note that 
the sequence starting mass transfer from a star at the terminal main-sequence defines a lower bound for 
the data points in Fig. H. It is difficult to see why CVs with somewhat evolved donors should dominate 
the population. One possibility is that unevolved systems are more susceptible to irradiation-driven mass 
transfer cycles (King et al. 1996), with a high state M too high to be recognisable as CVs. The scatter 
at short periods (3 — 5 hr) in turn seems to imply a wide range of secular mean M values, (10~ 9 — 10~ 8 
M yr _1 ). None of the orbital angular momentum loss mechanisms postulated so far generates this wide 
a range of M at a given P, while observational estimates for M always suggest that there is an even 
larger range (e.g. Patterson 1984, Warner 1995). 
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Fig. 5. Morphology of evolutionary tracks in the spectral type-orbital period diagram. Dashed: location of ZAMS 
secondaries. Solid curves: standard sequence (upper left panel), unevolved sequence with M — 1CP 7 Moyr" 1 
(upper right), evolved sequence (X c = 4 x 10~ 4 ) with M — 1.5 x 1CT 9 M Q yr _1 (lower left), and with M = 5 x 10~ 8 
Moyr - (lower right). 



If the M range reflects the range of the secular mean, then the appearance of a gap with relatively sharp 
edges is more difficult to understand than with a uniform evolution driven by angular momentum losses 
insensitive to system parameters (Kolb 1996; Stchle ct al. 1996). Yet the observed period distribution and 
its explanation by disrupted orbital braking cannot unambiguously rule out a large range of the secular 
mean transfer rate above the gap (Baraffe & Kolb 1999). The M range certainly has to level off rather 
steeply below 10 -9 M0yr _1 to avoid filling the gap with systems that hardly detach. On the other hand, 
systems with M^2x 10~ 9 M Q yr _1 undergo period bounce before they reach the upper edge of the gap, 
and before the donor becomes fully convective (at a mass M CO nv)- This bounce period increases with M, 
while Mconv decreases. Hence the detached phase of high M systems extends from periods longer than 
3 hr to periods shorter than 2 hr, and does not interfere with the "classical" gap. The time spent in 
the detached phase increases as well, in the case M = 1CP 8 M Q yr _1 it is longer than 8 Gyr (assuming 



angular momentum losses by gravitational radiation). Only a small fraction of such systems would reach 
the orbital period minimum at 80 min within the age of the Galaxy thus alleviating (but not solving) 
the problem of the missing period spike (see Kolb & Baraffe 1999) at this minimum period. 

Our conclusions from the mass-spectral type diagram and orbital period-spectral type diagram can 
be tested by a set of more accurate mass and radius determinations of CV secondaries. 

References 

Allard, F. 1999, in Very Low-Mass Stars and Brown Dwarfs in Stellar Clusters and Associations, Euroconference, 

La Palma 1998 
Allard, F., Hauschildt, P.H., Baraffe, I., Chabrier, G. 1996, APJ, 465, L123 
Allard F., Hauschildt P.H., Alexander D.R., Starrfield S. 1997, ARA&A, 35, 137 
Baraffe I., Chabrier G. 1996, ApJ, 461, L51 
Baraffe I., Kolb U. 1999, MNRAS, submitted 

Baraffe I., Chabrier G, Allard F., Hauschildt P.H. 1995, ApJ, 446, L35 
Baraffe I., Chabrier G, Allard F., Hauschildt P.H. 1997, A&A, 327, 1054 
Baraffe I., Chabrier G, Allard F., Hauschildt P.H. 1998, A&A, 337, 403 
Beuermann, K., Baraffe, I., Kolb, U., & Weichhold, M. 1998, A&A, 339, 518 
Chabrier G, Baraffe I. 1995, ApJ, 451, L29 
Chabrier G, Baraffe I. 1997, A&A, 327, 1039 
Kalogera V., Kolb U., King A.R. 1998, ApJ, 504, 967 
King A.R. 1988, QJRAS, 29, 1 
King A.R., Kolb U. 1997, ApJ, 481, 918 
King A.R., Frank J., Kolb U., Ritter H. 1996, ApJ, 467, 761 
Kolb U. 1996, in Cataclysmic Variables and Related Objects, ed. A. Evans, J.H. Wood, Dordrecht: Kluwer, IAU 

Coll. 158, 433 
Kolb U., Baraffe I. 1999, MNRAS, submitted 
Kolb U., Ritter H. 1992, A&A, 254, 213 
Kolb U., King A.R., Baraffe I., 1999, MNRAS, submitted 
de Kool M. 1992, A&A, 261, 188 
Mazzitelli I. 1989, ApJ, 340, 249 
Patterson J. 1984, ApJS, 54, 443 
Politano M. 1996, ApJ, 465, 338 
Ritter H., 1994, in Evolutionary Links in the Zoo of Interactive Binaries, ed. F. D'Antona, V. Caloi, C. Maceroni, 

F. Giovanelli, in Memorie della Societa Astronomica Italiana 65, 173 
Saumon, D., Chabrier, G, and VanHorn, H.M., 1995, ApJS, 99, 713 
Saumon, D., Chabrier, G, Wagner, D.J., Xie, X., 1998, Phys. Rev., submitted 
Smith D.A., Dhillon V.S. 1998, MNRAS, 301, 767 
Stehle R., Ritter H., Kolb U. 1996, MNRAS, 279, 581 

Warner B. 1995, Cataclysmic Variable Stars, Cambridge Astrophysics Series 28, Cambridge: CUP 
Whyte C, Eggleton P.P. 1980, MNRAS, 190, 801 



